Atomic collapse, Lorentz boosts, Klein scattering, and other quantum-relativistic 

phenomena in graphene 
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Electrons in graphene, behaving as massless relativistic Dirac particles, provide a new perspective 
on the relation between condensed matter and high-energy physics. We discuss atomic collapse, 
a novel state of superheavy atoms stripped of their discrete energy levels, which are transformed 
into resonant states. Charge impurities in graphene provide a convenient condensed matter system 
in which this effect can be explored. Relativistic dynamics also manifests itself in another system, 
graphene p-n junctions. We show how the transport problem in the presence of magnetic field can 
be solved with the help of a Lorentz transformation, and use it to investigate magnetotransport in 
p-n junctions. Finally, we review recent proposal to use Fabry-Perot resonances in p-n-p structures 
as a vehicle to investigate Klein scattering, another hallmark phenomenon of relativistic dynamics. 
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I. INTRODUCTION 

The unique properties of graphene, a single atomic 
layer of carbon [J], attract a lot of attention and inter- 
est from researchers from diverse fields. What makes 
this material special is its transport characteristics, such 
as high mobility and tunable carrier density, which of- 
fer new exciting opportunities for nanoelectronics . At 
the same time, charge carriers in graphene exhibit many 
unusual properties, posing interesting questions of fun- 
damental interest [3| . 

Undoped graphene is a semimetal in the solid state 
nomenclature, with conduction and valence bands joined 
together at the Fermi points, which are the symmetry 
points of the Brillouin zone known as K and K'. Near 
these points, the conduction band has linear dispersion, 
which can be modeled by the Dirac equation. Charge 
carriers in graphene are thus described as massless rela- 
tivistic fermions, with an effective "speed of light" v w 
lO^m/s, which is about 300 times less than c, the speed 
of light in vacuum. 

This description in terms of massless Dirac particles 
leads to an interesting analogy between the physics of 
graphene and that of high-energy relativistic particles, 
offering a possibility to observe quantum-relativistic phe- 
nomena in condensed matter experiments. Various inter- 
esting phenomena associated with relativistic fermions, 
such as the half-integer quantum Hall effect |4|(5j|, Klein 
tunneling HQ and scattering HQ, atomic collapse [lOl . 
HlllT^ . gauge fields and topological defects [isllT^ . can be 
translated directly into transport properties of graphene. 

Compared to other two-dimensional electron systems, 
the linear dispersion in graphene, e — ±f|p|, leads to 
much higher quantization energies. This is the case, 
for example, for the quantized energies of Landau lev- 
els. The Dirac Hamiltonian in a uniform magnetic field 
yields e„ = zLv\/2ehBri, n — 0, 1,2..., with an interlevel 
spacing that can exceed 1000 K in a field of about 6 Tesla. 
As a result, quantized Hall transport in graphene can be 



observed at room temperature [l5|. 

The "Diracness" of charge carriers plays a fundamental 
role in determining the electronic properties of graphene 
(Sj. In terms of the wavefunction amplitudes on the A 
and B sublattices of the carbon honeycomb lattice, tp = 
{^Pa,iPb), the Hamiltonian takes the form 
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p±=pi± ip2 (1) 



where C, = -l-l(— 1) for the point K{K'). Starting with 
the Schrodinger equation ihdtip = H'ip, and multiplying 
it by the matrix cr^ , we can write it in the canonical form 
of a relativistic Dirac equation for a free particle 



-/''df.ip = 0, = {vt, xi,X2)- 



(2) 



The 2x2 matrices 7^" satisfy the anticommutation rela- 
tions {7^,7"} = g^'^, where g^'^ is the metric tensor. 

The analogy with special relativity extends to prob- 
lems involving coupling of electrons to external electric 
or magnetic fields. Rewriting the Schrodinger equation 
ihdtip = [H{p — ea) -I- e(p] ip in the form ^ yields the 
Dirac equation in a background electromagnetic field, 

j'^ihd^ - ea^)ip ^ 0, = (i;"Viai,a2). (3) 

The potentials describe external fields E = — Vi^ and 
B ~ V X SL, pointing in-plane and out-of-plane, respec- 
tively, which can be either static or dynamical. The anal- 
ogy with 2-1-1 Quantum Electrodynamics (QED), man- 
ifest in Eq.([3]), helps to establish connections between 
various problems involving electrons in graphene with 
analogous problems in QED, described by the relativistic 
Dirac equation. 

One aspect of the relation with QED which is particu- 
larly intriguing is the large value of the dimensionless cou- 
pling strength, the analog of the fundamental fine struc- 
ture constant for electrons in graphene. Indeed, since 
c/v ~ 300, we have 
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where the dielectric constant k, describing the cumulative 
effect of screening by the substrate and intrinsic screening 
ingraphene, typically can take values between 3 and 10 
\yf\. The large value of the fine structure constant, a ^ 1, 
makes the QED-likc physics in graphene even more rich 
and interesting than in "natural" QED systems. 

As a cautionary remark we note that, while the anal- 
ogy with 2+1 QED is valid for coupling to arbitrary 
external fields, and is indeed useful in a variety of prob- 
lems, it does not extend to the dynamical EM field caused 
by charges in graphene. The basic reason for this differ- 
ence is the three-dimensional character of such EM fields, 
which originate on the charges in graphene, but spread 
out in the entire space outside the graphene plane. Hence 
the problem of electron interactions in graphene requires 
insights which are outside the realm of 2 + 1 QED [e.g., 
see Refs . [18111911201 and references therein]. 

In this article we review some of the recent work on 
graphene motivated by analogy with QED, focusing on 
the phenomena familiar from relativitic quantum me- 
chanics and on their connection to electronic transport. 

In Sec. |TT]we shall discuss the problem of atomic col- 
lapse, in which the large value of a is crucial. Atomic col- 
lapse is a phenomenon discovered theoretically in atomic 
physics, and predicted to occur in heavy atoms, hav- 
ing nuclear charge Z > 170. Because such large val- 
ues of Z are not realized in any of the known elements, 
this phenomenon has not been experimentally observed. 
As we shall see, in graphene the condition for collapse, 
Z > ~ 1, is much easier to realize. In particular, 
collapse can occur for the electronic states near c harg e 
impurities, with critical charge as low as = 2 [lollll| . 

In Sec. mil we consider another system for which the 
connection with relativistic dynamics turns out to be 
fruitful: graphene p-n junctions. We show that trans- 
port in a p-n junction in the presence of a magnetic field 
can be understood with the help of a Lorentz transfor- 
mation, which is a symmetry of the Dirac equation ([3]). 
This transformation, in which the velocity v replaces the 
speed of light, can be used to eliminate either the mag- 
netic field or the electric field, whichever is weaker. We 
compare this method to other approaches, and use it to 
find the field-dependent conductance of the system Q . 

Next, in Sec. llVi we discuss transport in p-n-p struc- 
tures. We focus on the ballistic regime, in which conduc- 
tance exhibits Fabry-Perot (FP) resonances due to inter- 
ference of electron waves refiected from two p-n bound- 
aries. This system is convenient for investigating the 
Klein scattering phenomenon, which is a manifestation 
of chirality conservation in the dynamics of relativistic 
Dirac particles. In particular, because the reflection coef- 
ficient vanishes for normally incident particles, the back- 
refiection amplitude, an analytic function of the incidence 
angle, changes sign at normal incidence. This sign change 
contributes a tt phase shift to the FP interference, which 
can be revealed by a half-a-period shift of FP fringes in- 
duced by a relatively weak magnetic field Q . These pre- 
dictions are in agreement with recent experiment '2l|. 



II. ATOMIC COLLAPSE AND 
SUPERCRITICAL CHARGE IMPURITIES 



As we mentioned above, graphene offers an opportu- 
nity to observe the behavior of relativistic particles in 
very strong fields, a regime which is very difficult to 
achieve in high-energy experiments because of the small 
value of the fine structure constant. This is exemplified 
by an interesting phenomenon, which will be discussed 
in this section, the collapse of a Dirac particle moving in 
a 1/r potential. 

To gain some intuition about this problem, let us start 
by recalling the explanation of stability of the hydro- 
gen atom in the framework of Bohr's quantum theory. 
Atomic stability results from quantum zero-point mo- 
tion of an electron which prevents it from falling on 
the nucleus. Indeed, starting from the nonrelativistic 
Schrodinger equation, and estimating the kinetic energy 
of an electron as i^kin = fi^ /2rap^ , where p is the charac- 
teristic radius of electron wavefunction, we can write the 
total energy as 



Z^ 
P 



(4) 



where Z is nuclear charge. As a function of p, the energy 
(PI) is dominated by the positive kinetic energy term at 
small p, and by the negative Coulomb term at large p, 
giving rise to a minimum at po = / {Ze^m), the Bohr's 
radius. The latter quantity determines the size of an 
atom. 

The effects of quantum zero-point motion which pre- 
vent an electron from falling onto the nucleus are less 
powerful in the relativistic regime. In this case, the bal- 
ance between kinetic and potential energy becomes more 
delicate, because the kinetic energy e = cp ^ ch/ p scales 
in the same way as the Coulomb energy. To see this 
more clearly, let us extend the reasoning used in (|4]) to 
a relativistic electron, described by the kinetic energy 
e{p) = {cpy + [mc^Y. Estimating p '-^ ?i/p, we again 
seak to minimize the sum of the kinetic and potential 
energies: 




(mc^y 



(5) 



Looking for extrema of this function, we find an equa- 
tion ^1-1- {mc/TiYp'^ — hc/Ze^, which has a solution 
only when Z < Tic/e^ « 137. At higher values of Z the 
effective potential ([U is a monotonic function varying be- 
tween E = mc^ at p = cx) and E = —oo at p = 0. Thus 
there is no stable solution for a hydrogenic atom at large 
nuclear charge. 

The difference between the high and low values of Z 
can be also seen directly in the exact solution of the Dirac 
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FIG. 1: a) Energy levels of superheavy atoms obtained from 
Dirac equation for Coulomb potential ~Ze^ /r, plotted as a 
function of = Za, where Z is nuclear charge, and a = /he 
is the fine structure constant. Energy is in the units of mc? . 
(b) Energy levels for Coulomb potential regularized on the 
nuclear radius. As Z increases, the discrete levels approach 
the continuum of negative-energy states and dive into it one 
by one at supercritical Z > 170 (from Ref.jU 



equation for a hydrogenic atom. 



En,j - 



-1/2 



vifith the quantum numbers n — 1,2. 



(6) 

= ±1,±2..., 



and a = e /Tic ~ 1/137. The expression under the square 
root becomes negative when Za > 1, rendering the eigen- 
values ([6]) with j — ±1 complex- valued. The unphysical 
complex energies indicate an intrinsic problem arising in 
the Dirac equation at supercritical Z > a~^. 

Mathematically speaking, this behavior can be under- 
stood as an effect of the 1/r singularity of the Coulomb 
potential, making the Dirac operator non-Hermitian and 
leading to the breakdown of the Dirac equation. One 
could replace the 1/r potential by a regularized poten- 
tial, rounded on a characteristic nuclear radius, which can 
be achieved by using a suitable nuclear form factor ^2^. 
However, it was found that, after the 1/r singularity is 
eliminated and the Dirac equation becomes well defined 
as a mathematical problem, the anomalous behavior at 
large Z persists. 

This problem was analyzed half a century ago 2^l23j , 
and it was predicted that in superheavy atoms with 
Z > 170 relativistic effects lead to a reconstruction of 
the Dirac vacuum. As the nuclear charge Z increases, it 
was found that the energies of discrete states approach 
the negative energy continuum, e < — mc^, and then dive 
into it, one after another (sec Fig[T]). After entering the 
continuum, discrete states turn into resonances with a 
finite lifetime, which can be described as resonant (or, 
quasistationary) states with complex energies. 

The finite lifetime of electronic states, which indicates 
an instability of superheavy atoms, puts a natural limit 
on the extent of the periodic table of elements. As such, 
the prediction of this phenomenon caused a lot of inter- 
est and excitement. Experiments on heavy-ion collisions 
probing this collapse have been suggested and attempted 



[2J|; however, the results were ambiguous and difficult to 
interpret. 

Graphene offers a completely new perspective on the 
problem of supercritical atoms, with charged impurities 
providing a natural realization of the Coulomb poten- 
tial. Charge carriers in graphene, which are described 
by a massless Dirac equation, mimic Dirac electrons in 
the ultrarelativistic regime, e 3> mc^. In this case, an 
analysis of the balance between the zero-point motion 
and Coulomb attraction similar to that leading to Eq. ^ 
yields 



hv Ze^ hv — Ze^ 

E{p) = = • 

P P P 



(7) 



For Z > a^^ — hv/e^ , the energy can be driven to arbi- 
trarily large negative values by letting the radius p tend 
to 0; thus we find a break-down or "collapse" of ultra- 
relativistic atoms in the supercritical regime Z > a"^. 
Due to the low value of the "speed of light" the critical 
charge in graphene is of the order of one [l3| • This makes 
charged impurities in graphene a very convenient system 
for experimental investigation of this phenomenon. 

Electronic states near a charged impurity are described 
by the two-dimensional Dirac-Kepler problem 
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Properties of the solutions of Eq.® depend on the di- 
mensionless parameter (3 — Ze^ /Uvk, exhibiting two dif- 
ferent regimes. For < 1/2, the behavior of states 
is consistent with scattering on a weak potential. For 
\(3\ > 1/2, there is an abrupt reconstruction of the Dirac 
vacuum around the 1 /r potential, which leads to the for- 
mation of resonant states lying in the continuum. These 
states can be understood as discrete electronic states, 
broadened into resonances by their coupling to the con- 
tinuum of hole states [ll|. Comparing this behavior 
to the predictions for superheavy atoms 2^12^ . we can 
identify these two regimes with 

no collapse {\f3\ < 1/2) and collapse {\f3\ > 1/2). (9) 

Note the critical value 1/2 instead of 1 foimd in three 
dimensions. 

The formation of resonances can be understood using 
scattering states of the Dirac-Kepler problem. Using po- 
lar coordinates, x + iy — pe"^, the general solution of 
Eq. ([8]) can be written as a sum of cylindrical waves with 
different angular quantum numbers [To| 



w{p)+v{p) 

[w{p) - v{p)y^ 



^l/2gi(m- 



(10) 

where s — ^ vn? — /J^, fc = —e/Tiv^ and m is a half- 
integer angular quantum number. Here w and v are 
the amplitudes of the incoming and outgoing cylindrical 
waves, satisfying a hypergeometric equation. Solutions 




Potential strength [i 

FIG. 2: (A) The dependence of the scattering phases S m on 
energy e = —hv/k, shown for several subcritical and super- 
critical values of /3 = Ze^ /hv (from Ref. [lol|). For subcritical 
/3, the phases Sm have no energy dependence (black line). For 
supercritical /3, the phases acquire energy dependence, which 
becomes stronger as |/3| increases. For /3 < — i (Coulomb 
attraction), the kinks in 5m{k), indicate formation of reso- 
nant states with complex energies, Ree < 0, corresponding to 
Breit-Wigner resonances in scattering. For /3 > | (Coulomb 
repulsion) , the dependence 5m vs. log k is approximately lin- 
ear, exhibiting no kinks. (B) The transport cross-section, 
Eg. ljlip . shown as a function of /3 for several energy values. 
Note the Fano-shaped resonances at negative /3 < — i, ap- 
pearing due to resonant states (kinks in (A)), and oscillations 
at /3 > i, resulting from the energy-dependent Sm- 
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FIG. 3: (a) Local density of states (|12|) calculated at a fixed 
distance p = 10"^ ro from the charged impurity, where ro is a 
short-distance parameter of the order of carbon lattice spacing 
(from Ref.(lQl]V Peaks in the LDOS, which appear at super- 
critical (3 and move to more negative energies at increasing 
j/3|, correspond to the resonant states, (b) Spatial map of 
the density of states, shown for several values of /3, with res- 
onances marked by white arrows (from Ref. [ill]). Note that 
the spatial width of the resonances decreases at they move 
to lower energies, Ap oc l/|e|, while the linewidth increases, 
7 oc \e\. The oscillatory structure at positive energies repre- 
sents standing waves with maxima at kp (n + |)7r, similar 
to those studied in carbon nanotubes [26l |. Energy is given in 
the units of £o = W~^Tiv/ro ~ 30 mV for ro = 0.2 nm. 



of this equation, combined with a boundary condition 
at short distances, can be used to construct scattering 
states. A "zigzag" boundary condition was used in 
with V'2(p) = at the radius p = ro on the order of 
carbon atom spacing. Other boundary conditions yield 
similar behavior. 

Scattering phases Sm{k), obtained from the asymptotic 
ratio -y/w — Q^^Sm{k)+2ikp ^j^^ ^ ^-^^ exhibit very different 
behavior in the two regimes illustrated in Fig[21 The 
resonant states with complex energies, formed at /3 < 
— i, appear as kinks in the dependence 6m vs. k. 

The resonant states should manifest themselves exper- 
imentally via various transport properties. In particular, 
the transport scattering cross-section 



Ctr = 1" 



sin^((5„i+i - dm), 



k ^ 

m=0 



(11) 



exhibits a resonant structure which is shown in Fig[2}3. 
In contrast, in the subcritical regime |/3| < ^ the differ- 
ences (Jm+i ~ <^m are energy- independent, which makes 
the cross-section scale as with a slight asymmetry in 
the prefactor for positive and negative /3 values, which 



was discussed in |25l |. 

The most direct signature of resonant states, however, 
can be seen in the local density of states (LDOS), 



v{e,p) 



Trhv 



(12) 



where tp{k,p) are the scattering states ([TU)) . This quan- 
tity can be directly measured by energy-resolved scan- 
ning tunneling spectroscopy probes. 

The energy dependence, Eq. p2)) . calculated at a fixed 
distance from the charged impurity, is illustrated in FiglJl 
Overall, the LDOS is an approximately linear function 
of energy, with oscillations at positive energies and sharp 
resonances at negative energies, which appear at super- 
critical p. The resonances, which occupy a spatial region 
of size p ~ hv/\e\ near a charged impurity, provide a 
clear signature of atomic collapse accessible with scan- 
ning probe techniques. 

To estimate the critical value of impurity charge Zc — 
l/2a, we need to find a = e'^/hvK taking into account 
screening. Screening due to the dielectric substrate with 
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dielectric constant e is described byK=(e+l)/2. More 
interesting, however, is the intrinsic screening due to po- 
larization of the graphene Fermi sea. Using random- 
phase approximation (RPA) for the dielectric function 
of undoped graphene, krpa ~ 5, and ignoring the effect 
of dielectric substrate, we find that Z = 1 is subcritical, 
while Z = 2 is supercritical llOll . This conclusion is in 
agreement with recent work [27| which analyzes screen- 
ing of charged impurity as polarization of Dirac vacuum 
in a strong Coulomb field (see also Refs. 20l[28t ). 

It was noted in Ref. [10| that linear screening breaks 
down in the supercritical case, for undoped graphene. 
The reason is that the vacuum polarization in the pres- 
ence of supercritcal l/r potential, found from the solu- 
tion of the noninteracting problem which we considered 
above, has a profile. This leads to a log divergence 
of the polarization charge, and thus to overscreening, un- 
less the effects of interaction within the polarization cloud 
are taken into account. It was shown with the help of a 
renormalization group argument that the true value of 
polarization charge is —{Z — Zc), i.e. the polarization 
cloud compensates the excess part of the impurity charge 
and brings it down to the critical value. 

The profile of the induced polarization can be found 
approximately for e^/hv <C 1 [l^, and also for e^/hv ^ 1 
[29|. The effective potential, resulting from adding a con- 
tribution of polarization charge to the bare 1 /r potential, 
deviates from Coulomb potential, and thus the above so- 
lution of the Dirac-Kepler problem strictly speaking is 
not applicable. Still, the qualitative picture described 
above, with resonant states appearing in the Dirac con- 
tinuum at Z > Zc remains valid even in the nonlinear 
screening regime. 



III. LORENTZ TRANSFORMATION AND 
COLLIMATED TRANSMISSION IN P-N 
JUNCTIONS. 

A basic property of massless Dirac particles is conser- 
vation of chirality, defined as the projection of pseudospin 
on the velocity vector. Conserved chirality leads to unim- 
peded tunneling through arbitrarily high barriers (the so- 
called Klein phenomenon) [s^l; particles with conserved 
chirality, incident normally on a barrier, exhibit perfect 
transmission and zero reflection [^Q . The Klein effect al- 
lows electrons to leak out of any trap fssiiBQi] , and makes 
electrostatic confinement in graphene more difficult than, 
e.g., in gated semiconductor systems. Klein transmission 
is also an essential part of theoretical understandingof 
transport in systems such as graphene p-n junctions [3|. 

The range of angles in which transmission can occur 
depends on the microscopic details of a potential bar- 
rier 0. For smooth barriers, such as those in which the 
potential varies on a scale large compared to the Fermi 
wavelength, perfect transmission is restricted to a narrow 
range of angles. This is the case for p-n junctions created 
by electrostatic gates, which are typically placed few tens 



of nanometers above graphene plane [3 

Collimated transmission in such systems can be ana- 
lyzed by replacing the externally imposed electrostatic 
field in the region where Klein tunneling occurs by a uni- 
form in-plane electric field 7]. For massless Dirac par- 
ticles moving in a uniform in-plane electric field -E |j x, 
the transmission coefficient, as a function of electron mo- 
mentum projection on the p-n interface, is given by 



t{Py 



X = 7rv/h\eE\ 



(13) 



The range of py in which transmission can occur becomes 
very narrow for small i.e. for smooth potentials. 
What happens to the collimated Klein transmission 
in the presence of a magnetic field? The Dirac 
equation ([3]), describing an electron in graphene in the 
presence of external fields, exhibits two qualitatively dif- 
ferent regimes, depending on the relative strength of the 
E and B fields. For weak magnetic fields, \B\ < \E\/v, 
the picture of collimated transmission described in Ref. [3| 
remains essentially unchanged, with perfect transmission 
occurring for the incidence angle which is a function of 
magnetic field. 



00 = arcsin(t>B/£;), ti9o) = 1. 



(14) 



In contrast, at stronger magnetic field, v\B\ > \E\, elec- 
tronic states are described by quantized Landau levels 
with position-dependent energies [30]. In this regime, in 
a clean system, transport can occur only perpendicular 
to the electric field. 

There are several ways to understand the origin of the 
two regimes. Perhaps most easily and elegantly, it can 
be seen by employing Lorentz invariance of the equa- 
tions of motion, Eq.([3]), with respect to Lorentz trans- 
formations in which the graphene dispersion velocity v 
plays the role of the speed of light. There is a basic 
fact, established in special relativity, that for perpendic- 
ular electric and magnetic fields, a Lorentz boost can 
be performed such that in the moving frame one of the 
fields vanishes, while the other field remains nonzero, ex- 
periencing Lorentz contraction. The field which can be 
eliminated by such Lorentz boost is the weaker of the 
two. Thus, for v\B\ < \E\, one can eliminate the mag- 
netic field by performing a boost to a moving frame with 
velocity u — v^B/E. For v\B\ > \E\, the electric field 
can be eliminated by a boost with velocity u — E / B. In 
the first case, with B zero and E nonzero in the moving 
frame, we have a situation identical to that studied in 
Ref. Q. In the second case, with i? = in the moving 
frame, and B finite, electronic states are derived from 
quantized Landau levels. 

It may be useful to point out that the situation with 
nonrelativistic, Galilean transformations is quite differ- 
ent. In the presence of perpendicular electric and mag- 
netic fields, a Galilean transformation can be used to 
eliminate E by performing a boost in the moving frame 
with velocity u = E/B. It is not possible, however, to 
eliminate magnetic field by a Galilean transformation. 
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Given this property of Gahlean transformations, one 
may be suspicious about drawing conclusions from 
pseudo-Lorentz transformations, which appear to be a 
kind of mathematical trick, rather than a genuine sym- 
metry of the system. To relieve this concern, we now 
present another argument leading to the two regimes, 
which does not rely on Lorentz invariance (following a 
suggestion by F. D. M. Haldane). 

Let us consider classical trajectories of a charged parti- 
cle with kinetic energy e(p), which is moving in the fields 
E and B. This problem can be described by canonical 
equations of motion with a Hamiltonian 

Hipjr) ^ e{p) - eEx, p = p - eA, A= {0,Bx), 

where we have chosen the gauge so that the vector po- 
tential is parallel to the vector Vd = E x B/i?^, which 
defines classical drift velocity. With this choice of A, and 
taking into account conservation of the component of mo- 
mentum p2 perpendicular to Vd , it is straightforward to 
write an equation for particle trajectory in momentum 
space, given by the energy integral 



e(p) - Vd.p 



where we used the relation P2 = P2 — eBx to express x 
through p2, absorbing constant p2 in Eq. In the case of 
massless Dirac particles, with the kinetic energy £(p) — 
±w|p|, we can write the trajectory in polar coordinates 
in momentum space as 



p{e) = 



£0 



ztv — Ud cos ( 



where 9 is the angle measured relative to Vd. This for- 
mula is nothing else than an equation for conical sections, 
familiar from the theory of planetary motion (an electron 
is orbiting the Dirac point in a manner similar to a planet, 
or a comet, orbiting the Sun). Depending on the relative 
magnitude of v and v^, there are two cases, 

Vd = E / B > V (hyperbola), Vd = E/ B < v (ellipse), 

which coincide with the two regimes identified using 
Lorentz transformations. The trajectories are open for 
Vd > V, and closed for Vd < v, corresponding, respec- 
tively, to motion that originates at infinity and ends at 
infinity, and to periodic cyclotron motion. 

Now we proceed to analyze transmission in the regime 
v\B\ < \E\. As discussed above, we can eliminate i? by a 
Lorentz boost with rapidity /? = —vB/E parallel to Vd: 



7 7/3 
1 
7/3 7 



1 



7 = 



(15) 



In the new frame we have B' = 0, E' = E/j. 

Since B' — 0, we can use Eq. fT^ . with E replaced 
by giving the transmission coefhcient as a function 



of momentum Py parallel to Vd- Expressing Py and E' 
through the quantities in the lab frame, we obtain 



TiPy) 



e/v. 



(16) 



Introducing the incidence angle, vpy = £sin0, we find 
the value 6 = aicsm(vB / E) for which the transmission 
is perfect, in agreement with Eq. (|14p . The angular size 
of transmitted beam is a function of B, such that the 
beam becomes more coUimated at higher B. 

In passing between the moving and lab frames we used 
the fact that the transmission coefficient, Eq. p^ . is a 
scalar with respect to Lorentz transformations, Eg. psp . 
This is indeed true because transmission and reflection 
at the p-n interface is interpreted in the same way by all 
observers moving with velocity u || Vd. 

The conductance can be found by integrating the 
transmission (|16p over py , which gives Q 



GiB) 



Go 




[i-ByBi 



2^3/4 



B<Bo 
B>Bo 



(17) 



where Bq — \E\/v = Tr/h\e\X and Go is the conductance 
of a ballistic p-n junction at S = 0, proportional to the 
length of the p-n interface, which was found in 0. 

The critical field value Bq can be estimated from den- 
sity profile in the p-n junction region, which can be found 
from electrostatic simulation, as discussed in Refs.fUlH. 
Using the relation n = kp/n, we estimate the electric 
field in the junction as eE w 2hvkp/d = 2hvy/Trn/d, 
where d is the width of the density step. This gives 
Bq « '^QyJripK / d, where $0 = h/e is the flux quan- 
tum. For n = 10^^ cm"^ and d = 100 nm we estimate 
Bo ~ 2.3 T. 

Caution should be used when the result p7|) . describ- 
ing the field dependence of conductance in the ballistic 
regime, is applied to real systems, such as those stud- 
ied in Refs. |3lll3^ . First, our derivation only accounts 
for the current flowing through the entire cross-section 
of the system, and not for, e.g., the edge currents. The 
latter become important at strong fields, B > Bo, when 
the system enters the quantized Hall regime, giving a 
contribution to conductance of the order of e^ /h. 

In addition, in studying a p-n junction in the presence 
of disorder, one should account for net contribution of 
disordered regions to total resistance. A series resistance 
model was proposed in Ref. [33l | to describe this situation. 
However, the magnetic field dependence in the disorder- 
dominated transport regime has not yet been analyzed. 



IV. KLEIN SCATTERING AND FABRY-PEROT 
RESONANCES IN P-N-P STRUCTURES 

Theoretically, Klein scattering is regarded as a funda- 
mental mechanism of transport through potential barri- 
ers and p-n junctions in graphene 00]. Until recently, 
however, this theoretical picture was disconnected from 
experiment. The manifestations of Klein scattering that 
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FIG. 4: Transmission coefficient of p-n-p structure, obtained 
from numerical solution of the Dirac equation with potential 
U{x) = ax^ — £, plotted as a function of the component of 
electron momentum parallel to the p-n boundaries Py and 
potential depth e. At zero magnetic field (a), transmission 
exhibits fringes with a phase which is nearly independent of 
Py. At finite magnetic field (b), fringe contrast reverses its 
sign on the parabola (black line) , which marks the boundaries 
of the interval (l22j (from Ref. 9]). Here e* = [ah'^v^y/^, 
Pt = e,/v. 



have been discussed so far, such as collimated transmis- 
sion which becomes perfect at normal incidence, can be 
easily masked by scattering on disorder. Because disorder 
is an integral part of any realistic system, these predic- 
tions may be difficult to verify using existing transport 
data. 

Recently, however, another manifestation of the Klein 
phenomenon in transport was discussed, which may of- 
fer a more direct experimental signature Q. The idea 
is to refocus attention from Klein transmission to Klein 
backreflection, which must vanish at normal incidence. 
Since the reflection amplitude is a smooth function of in- 
cidence angle 0, which vanishes at = 0, it must change 
sign when 9 varies from positive to negative values. Fur- 
thermore, because the sign change translates into a tt 
phase shift of the reflection phase, it can be revealed by 
Fabry-Perot interference in a p-n-p structure. 

The simplest regime to analyze is ballistic transport 
through a p-n-p structure. Transmission through two 
parallel p-n interfaces, which we label 1 and 2, can be 
described by the Fabry-Perot model, 



TiPy) = 



1^2 



1 1 - y/nr^e 



iA9 



(18) 



Weak B: 

zero not enclosed, AG = 




Strong B: 

zero enclosed, A9 = ji 
py i 



px 



FIG. 5: Berry phase for periodic orbit of an electron bounc- 
ing between p-n boundaries. As B increases, with py kept 
constant, the trajectory evolves from a skinny oval to a less 
skinny oval. For typical py the trajectory does not enclose 
the origin at small B, giving zero Berry phase. At higher B, 
when the origin is enclosed, the Berry phase equals tt. 



where ^1(2) and ri(2) = 1 



(2) 



^1(2) a-re the transmission and 



reflection coefficients for each interface. The phase AO is 
a sum of the WKB phase and the phases of the reflection 
amplitudes. 



Ae= - Px{x')dx' + A9i + A02 



(19) 



Here A9i(2) are the backreflection phases for the inter- 
faces 1 and 2, exhibiting a 7r-jump at zero incidence angle. 

One could argue that, since the incidence angles for 
two parallel p-n interfaces are the same, the two 7r-jumps 
cancel each other in the net phase A61 + A92 , making the 
interference in insensitive to this effect. However, as 
pointed out in Ref. , this cancellation can be eliminated 
by applying a magnetic field, which curves electron tra- 
jectories and makes the incidence angles unequal. 

The effect of the magnetic field can be taken into ac- 
count by using the vector potential A — (0, Bx). In this 
gauge, the y-component of electron canonical momen- 
tum (i.e. the component parallel to the p-n interfaces) is 
conserved, and hence one can write the y-component of 
electron kinetic momentum as 



Py{x) ^ Py - eBx, 



(20) 



where the canonical momentum py is a constant of mo- 
tion that labels different scattering states. Using this 
relation, and evaluating the transmission coefhcients on 
each interface x — Xi(^2) with the help of Eq. (jl3p . we 
find ii(2) — e^'^P!'(^i(2)), From this, we obtain reflection 
amplitudes 



sgn ^(2^1(2))] e''--^Py'> y'l_e-^P?(-i(2)), (21) 

where we factored out the sign, responsible for the phase 
jump, and a regular part of the phase e*^'^"*, as follows 
from analyticity in py. 

In the model used in Ref.fol], the gate potential was 
described by a parabola U (x) = ax^ — e. In this case, e > 
creates p-n interfaces at a;i(2) = i-y/e/a. Linearizing 
U{x) near a;i(2), we find the parameter A = ^(ae)"^/^. 
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Energy s/s, 

FIG. 6: Resistance of a p-n-p structure, obtained from trans- 
mission, shown in FigU FP interference fringes in resistance 
as a function of potential depth e exhibit half-a-period shift 
when magnetic field is applied. The units for magnetic field 
are B. = {h/e)(a/nvf'^ (from Ref.Q). 



The sign jumps in expressions (|2ip occur when 
Py(xi{2)) = Oj where semiclassical electron trajectories 
approach the boundaries 1 and 2 at normal incidence. 
We can thus identify a range of values 



eBxi < py < eBx2 



(22) 



for which the reflection amplitudes have opposite signs. 
For these values of Pj,, the FP interference has a tt phase 
shift compared to what could be inferred from the WKB 
phase. 

Additional insight into the origin of the tt phase shift 
induced by magnetic field in FP interference is provided 
by a Berry phase argument, which focuses on the proper- 
ties of quasiclassical trajectories. There is an analogy be- 
tween the Dirac Hamiltonian ([T]) and the Zeeman Hamil- 
tonian of a spin 1/2 in a time varying magnetic field, 
with the pseudospin, associated with A and B sublat- 
tices, playing the role of spin. If particle motion is pe- 
riodic, over each period the wave function gains a phase 
given by a half of the solid angle swept by the effective 
magnetic field bcft(i) = wp(i)- In the case of cyclotron 
motion, this phase, combined together with the WKB 
phase, was used to explain half-integer quantization of 
Landau levels 

Here we apply a similar argument to the trajectory of 
an electron bouncing between p-n boundaries. It is in- 
structive to do the analysis in momentum space, focusing 



on the evolution with increasing B. At small B, since Py 
is nearly constant (see Eq. ipO)) ). the trajectory maps out 
a skinny oval, shown in Fig|5l For typical values of Py 
this oval does not enclose the origin, and thus the Berry 
phase associated with it is zero. For larger B the oval 
height increases, until eventually it encloses the origin. 
At this point the Berry's phase becomes equal to tt. 

This behavior can also be clearly seen in FigHl in 
which the transmission coefficient, obtained from nu- 
merical solution of the Dirac equation with potential 
U{x) = ax"^ — e, is displayed. Transmission exhibits FP 
interference fringes in the e direction both at zero field 
(panel a) and at finite field (panel b). In the latter case 
the fringe contrast also exhibits sign reversal across the 
black line. This line is the parabola py = ^eB^feJa^ 
which marks the boundaries of the interval (|22|) , and sep- 
arates regions with opposite phase contrast. 

Contrast reversal directly manifests itself in conduc- 
tance, which can be found by integrating transmission 
over Py. Since the fringes in transmission are nearly 
vertical (Fig[4]), the integral over py yields an oscilla- 
tory function of e. As the B field increases, making 
the region with inverted contrast wider, the oscillations 
in conductance shift by approximately one half-period 
(see FigE]). The shift occurs in relatively weak fields 
B - = {TL/e){a/TLvfl^, which gives a few hundred of 
milliTesla for realistic parameter values. 

Fabry-Perot oscillations in a ballistic p-n-p structure 
were reported in recent experiment [2l| . The behavior 
of the FP fringe contrast under applied magnetic field is 
reminiscent of that in Figl6l with a shift of approximately 
half-a-period observed in the fields of about 0.2 — 0.5 T. 
The data [2l[ also exhibits an interesting crossover be- 
tween FP resonances at i? = to Shubnikov-deHaas os- 
cillations at high i?, which indicates that momentum- 
conserving transport through the structure, which is 
responsible for FP resonances, gives way to disorder- 
dominated transport at high magnetic fields. 

In summary, graphene is a fertile system in which many 
interesting quantum-relativistic phenomena, or their con- 
densed matter analogs, can be explored. As a material 
amenable to various experimental techniques, graphene 
offers new interesting connections between condensed 
matter and high-energy physics. 

We thank F. D. M. Haldane, P. Kim, A. K. Savchenko, 
and A. Young for useful discussions. 
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